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Lattices in three dimensions are oft studied from the "reciprocal space" perspective of diffraction. 
Today, the full lattice of a crystal can often be inferred from direct-space information about three 
sets of non-paraUel lattice planes. Such data can come from electron-phase or Z contrast images 
taken at two tilts, provided that one image shows two non-parallel lattice periodicities, and the other 
shows a periodicity not coplanar with the first two. We outline here protocols for measuring the 
3D parameters of cubic lattice types in this way. For randomly-oriented nanocrystals with cell side 
greater than twice the continuous transfer limit, orthogonal ±15° and ±10° tilt ranges might allow 
one to measure 3D parameters of all such lattice types in a specimen from only two well-chosen 
images. The strategy is illustrated by measuring the lattice parameters of a 10 nm WCi-i crystal 
in a plasma-enhanced chemical- vapor deposited thin film. 



I. INTRODUCTION 



II. CALCULATIONS 



If 10 nm crystals presented themselves to our percep- 
tions as 10 cm hand specimens, new rules for direct-space 
crystallography might have emerged. For example, if by 
tilting the specimen one could locate a set of 4-fold cross- 
fringes, from which a tilt (at 45° to those fringes) by 
35.3° just brings one to a new set of fringes 15.5 percent 
larger in spacing, then the hand specimen is likely face- 
centered-cubic. It's reciprocal lattice definitely includes a 
body-centered-cubic array like that characteristic of face 
centered-cubic-crystals. 



Here we discuss the geometry of crystals from the per- 
spective of lattice imaging in direct space. The required 
instrumentation consists of a TEM able to deliver phase 
or Z contrast lattice images of desired periodicities (e.g. 
spacings down to half the unit cell side for cubic crys- 
tals), and a specimen stage with adequate tilt (e.g. two- 
axes with a combined tilt range of ±18°). For crystals 
with lattice spacings of 0.25nm and larger, many analyt- 
ical TEM's will work, while a high resolution microscope 
with continuous contrast transfer to spatial frequencies 
beyond 1/(0. 2nm) can do this for most crystals. We 
demonstrate the process experimentally by determining 
the lattice parameters of a tungsten carbide nano-crystal 
using a Philips EM430ST TEM. Appropriately orienting 
the crystal, so as to reveal its three-dimensional struc- 
ture in images, is a key part of the experimental de- 
sign and will be discussed in detail. In the process, 
strategies for supporting on-line electron-crystallography, 
for three-dimensional lattice-correlation darkfield studies 
of nanocrystalline and paracrystalline materials, and for 
stereo-diffraction analyses, are suggested as well. 
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A. Overview of a "tilt protocol" in action 

For the stereo lattice-imaging strategy discussed here, 
low Miller index (hence large) lattice spacings are both 
easier to see, and more diagnostic of the lattice. Ori- 
entation changes directed toward the detection of such 
spacings are needed. In this section, tilt protocols op- 
timized for getting 3D data from one abundant class of 
lattice types (namely cubic crystals) are surveyed. We 
begin with a list of candidate lattice types (e.g. f.c.c. 
or h.c.p.) based on prior compositional, diffraction, or 
imaging data. 

Three non-coplanar reciprocal lattice vectors seen 
along 2 different zone axes are sufficient for inferring a 
subset of the 3D reciprocal lattice of a single crystal. Of- 
ten these are adequate to infer the whole reciprocal lat- 
tice. Hence the goal of our experimental design is to look 
for at least 3 non-coplanar lattice spatial frequencies, in 
two or more images. We prefer images with "aberration 
limits" smaller than the analyzed spacings, to lessen 
chances of missing other comparable (or larger) spacings 
in the exit-surface wavefield. So as to tilt from one zone 
to another, the crystal must also be oriented so that the 
desired beam orientations are accessible within the tilt 
limits of the microscope. 

To illustrate, we exploit the fact (considered more fully 
below) that all cubic crystals will provide data on their 
three dimensional lattice parameters if imaged down se- 
lected zones separated by 35.3°, provided that spacings 
at least as large as half the cell side a are reported in 
the images. With the images discussed here we expect to 
"cast a net" for 3D data on any cubic crystals whose cell 
side a is larger than 2 x 0.19nm — 0.38rim. More than 
85% of the cubic close packed crystals and nearly 40% of 
the elemental b.c.c. crystals tabulated in Wyckoff for 
example, meet this requirement, as of course would most 
cubic crystals with asymmetric units comprised of more 
than one atom. 

Although 35.3° is too far for the eucentric tilt axis in 
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our microscope, combining two tilts gives us a range of 
35.6°. Two images therefore can be taken at orientations 
35.3° apart, namely at {-ffi = 15.0°, 1^2 = 9.7°) and = 
— 15.0°, 1)2 = —9.7°), respectively, where di and i?2 are 
goniometer readings on our double tilt holder. This yields 
an "effective" tilt axis that runs perpendicular to the 
electron beam. Its azimuth is 123.5° in the xy plane of 
our images. The coordinate system used will be discussed 
in more detail later. For the special case of fee crystals, 
the zones of nterest are [001] and [112]. Since the (220) 
lattice planes are parallel to both desired zones, the tilt 
must be along these planes. That is, (200) fringes seen 
down a 4- fold symmetric [001] zone must make an angle 
approaching 45° with respect to the effective tilt axis. 
Nearly a third of the randomly-oriented crystals showing 
[001]-zone fringes will be sufficiently close 

The experimental results were unambiguous. We found 
many four-fold symmetric images having spacings con- 
sistent with WCi„2:- This has an f.c.c. lattice with 
a — 0.4248nm j|, When such a zone was found with 
fringes making 45° to the effective tilt axis in the first 
image, a new spacing was seen in the second image mak- 
ing a 3D lattice parameter measurement possible. The 
experiment is illustrated in Figure ^ 



B. Experimental designs 

Here we seek three non-coplanar periodicities from two 
images, although the analysis also works if they are dis- 
covered singly in three images. Given Miller indices 
(hikili) and (/i2A:2?2) for any two periodicities (i.e. vec- 
tors gi and g2, respectively, in the reciprocal lattice) of 
any crystal, first find zone indices [uavawa] of the beam 
direction = gi x g2 needed to view both spacings in 
one image. The axis for the smallest tilt that will make 
the beam orthogonal to a third periodicity with indices 
(h^k^lg), and reciprocal lattice vector g3, may then be 
defined by the vector Vt = gs x r^. Lastly, zone indices 
[ubVbWb] for the beam after the specimen has been tilted 
around this axis so as to image the third periodicity, may 
be obtained from the expression = vt x ga. Note here 
that we treat the Bragg angle for electrons as small (i.e. 
less than one degree). Thus the actual tilt required will 
be a fraction of a degree less. 

Although these cross product calculations can be done 
by first converting for example to "c-axis" cartesian co- 
ordinates 1^, the simplest determination of needed pa- 
rameters is perhaps done using the metric matrix G of a 
prospective lattice [||: 



G: 



a»a a»b a»c 
b»a b»b b»c 
c»a c»b c»c 



(1) 



If we denote row vectors formed from Miller (or lattice) 
indices as {ijk\, and column vectors as \ijk), then the 
zone A indices obey gi • = {hikili\uAVAWA) = 



and g2 • = {h2k2l2\uAVAWA) = 0. From these two 
equations, [uy^Wy^WA] follows except for a multiplicative 
constant which is not important. Similarly, the (possi- 
bly irrational) Miller indices of the tilt axis {ht,kt,lt) 
may be determined from gt • r^i = {htktlt\uAVAWA) ~ 
and vt • g3 = = {htktlt\G-^\hzkzh) §. Only in this 
fourth equality does G affect the calculation, and for 
cubic crystals it then simply offers a multiplying con- 
stant. Finally, the zone B indices follow (to within a 
factor) simply from gt • = {htktlt\uBVBWB) =^ and 
g3 • rs = (hskshluBVBWB) = 0. 

Two parameters which determine the attractiveness 
and feasibility of a given experiment are the spatial res- 
olution, and range of specimen tilts, that the microscope 
is able to provide. For a given lattice type, it is useful 
to: (i) go through the hst of all pairs of periodicities cal- 
culating the tilt between the zone associated with that 
pair and any third spacing of possible interest, and (ii) 
rank the findings according to the minimum-spacing that 
must be resolved, and the range-of-tilt that the specimen 
undergoes. This has been done for face-centered, body- 
centered, and simple-cubic lattices, and the results illus- 
trated in Figure]^, with examples for fee and bcc spacings 
no less than half of the cell side illustrated in Figure |^. 
Note that this calculation needs to be done only once for 
each unit cell shape. Factors like the multiplicity of a 
given zone type might also figure into the design of ex- 
periments with randomly-oriented crystals, although we 
have not considered them here. 

High tilts can be used to lower measured spacing un- 
certainties, in directions perpendicular to the electron 
microscope specimen plane [||. Hence the protocols of 
interest for a given experiment may be those which ap- 
proach goniometer tilt limits, or at least the limits of 
specimen tiltability. 

Concerning the resolution limit to use, we suggest the 
spacing associated with the end of the first transfer func- 
tion passband in the micrograph of interest, sometimes 
inferrable from regions in the image showing disordered 
material. Even in this case, possible thickness and mis- 
orientation effects warrant caution |||, |l^ . 

One may of course explore spatial frequencies in the 
specimen up to the microscope information limit. How- 
ever, spherical aberration zeros in the transfer function 
introduce the possibility that the microscope will sup- 
press some spatial frequencies present in the subspecimen 
electron wavefield. To lessen this problem, HREM im- 
ages taken at different focus settings could be compared, 
if the foci were chosen so that spatial frequencies lost at 
one defocus are likely recorded at another. This would al- 
low measurement of three non-coplanar reciprocal lattice 
vectors, without missing any whose (reciprocal) length is 
shorter than the longest among those three. 

Lastly, of course, the protocol chosen may depend also 
on the specimen. For an image field containing hundreds 
of non-overlapping but randomly-oriented nano-crystals, 
only two micrographs could allow one to measure the 
three-dimensional lattice parameters of all cubic crys- 
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FIG. 1: Tilt protocols for determining the lattice parame- 
ters of face-centered, body-centered, and simple cubic crystals 
from two high resolution images, plotted according to the re- 
quired tilt range, and required resolution limit in units of the 
cubic cell side. Superimposed on this plot are the tilt range, 
and resolution in units of the unit cell side a = 0.4248nm of 
WCi-x, required for the experiment reported here. 



tal types present with cell sides a greater than 2ra. On 
the other hand, for a single crystal specimen of unknown 
structure, both a great deal of tilt range, and consider- 
able trial and error tilting (or guess work based on lattice 
models) might be required before a single set of three in- 
dexable non-coplanar spacings is found. 



C. Inferring 3D reciprocal lattice vectors from 
micrographs 

Consider a specimen stage with two orthogonal tilt 
axes (and associated rotation matrices) Ti and T2, both 
perpendicular to the electron beam (the second only so 
when the first axis is set at zero tilt). When the specimen 
is untilted {-di = 0°, ^2 — 0°), vectors in the recipro- 
cal lattice of the specimen may be described, in coordi- 
nates referenced to the microscope, as column vectors |g) . 
When this reciprocal lattice vector is tilted to intersect 
the Ewald sphere by some double tilt in the sequence 
T2{'&2) then ri(i9i), it's presence may be inferred from 
diffraction patterns or micrograph power spectra. In our 
fixed coordinate system, \g) has become \gm), where the 
m means that \gm) is associated with the lattice period- 
icity dm = 1/ g-m recorded on a micrograph. Using matrix 
notation, we might then write: 



|g™) =ri(i?i)T2 (1^2)1 g) 



(2) 





Components of gm may be determined from the polar 
coordinates (5, Lp) of a spot in the power spectrum of a 
recorded image, following: 



(3) 



where g is the length of the diffraction vector (e.g. in 
reciprocal nm) and Lp is its azimuth corrected for lens 
rotation. 

Hence we can calculate the "untilted-coordinates" 
of reciprocal lattice objects at inferred experimen- 
tally from micrographs, using 

\g) ^T-\d2)T-\d,)\^m) ^ A{d^,d2)\^m), (4) 
where we've defined: 



A{dud2) = T~\^2)T^\^i 



(5) 



The resulting xyz coordinates of reciprocal lattice fea- 
tures g, associated with the crystal while in the untilted 
goniometer specimen orientation, but referenceable from 
micrographs taken at any orientation, provide the lan- 
guage we use for speaking of our measurements in three 
dimensions. 



D. Calculating lattice parameters 

Given 3D cartesian coordinates of "points" in the re- 
ciprocal lattice of a crystal, we are in much the same situ- 
ation as if we had diffraction patterns of the crystal from 
two directions containing three (or more) non-coplanar 
spots. Hence methods for stereo-analysis of diffraction 
data [|[ |ri| can be used at this point. We summarize 
in this context briefly. Given measured reciprocal lattice 
vector coordinates, a natural next step is to infer a basis 
triplet for the crystal's reciprocal lattice. Three alter- 
nate paths to this basis triplet might be referred to as 
"matching", "building" [0, and "presumed" ||]. 

Given an experimental basis triplet from any of these 
sources, lattice parameters (a, 5, c, a, /3, 7), goniometer 
settings for other zones, and many other things follow 
simply from the oriented triplet matrix defined below: 



W = 



ax 


ay 


az 






b* 


c* ' 

^x 




by 






a*y 


K 


C* 




Cy 


Cz 




_ a-l 


K 


C* 



(6) 



Given W, for example, cartesian coordinates in the 
microscope for any direct lattice vector with indices 
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[uuTii] follow from |r) = W\uvw), while cartesian co- 
ordinates for the reciprocal lattice vector with indices 
(hkl) may be predicted from {g\ ~ {hkl\W^^. These 
rules of course include instructions for calculating basis 
vectors of the lattice, such as a = [100], and recipro- 
cal lattice, such as b* = (010), and the angles between. 
Moreover, the oriented cartesian triplet W is simply re- 
lated to the metric matrix for the lattice in equation 
|l| by G = WW^. From G, of course, all the familiar 
orientation-independent properties of the lattice follow 
0, including cell volume Vceii — \/\G\, Miller/lattice 
vector dot products {ghki\i'uvw) = {hkl\uvw), recipro- 
cal lattice vector and interplanar spacing magnitudes 
9hki = '^/'^hki — {hkl\G^^\hkl), lattice vector magnitudes 
''^uvw — {uvw\G\uvw) , reciprocal lattice dot products 
(51I52) = {hikili\G~'^\h2k2l2) , interspot angles 6*12 = 
cos"^ [(gi|ff2)/(5iff2)], etc. 

Even before a basis triplet is selected, indexing of ob- 
served reciprocal lattice vectors g can be attempted by 
matching spacings and interspot angles to candidate lat- 
tices. Because of the uniqueness of non-coplanar triplets 
in three dimensions (providing at least a significant sub- 
set of the whole reciprocal lattice) , the matches are very 
discriminating (even for low-symmetry lattices) relative 
to similar analyses from 2D data, i.e. from a only a sin- 
gle image or diffraction pattern. After a basis triplet is 
selected, the indices (hkl) of any observed spot {g\ in a 
diffraction pattern follows from {hkl\ ~ {g\W. Similarly, 
indices [uvw] of any observed lattice periodicity |r) in an 
image follow from \uvw) — W~^\r), once a basis triplet 
is in hand. Given the triplet one can similarly calculate 
goniometer settings to align the beam with any other 
crystallographic zone of interest. 

III. THE EXPERIMENTAL SETUP 

A. Instruments 

The Phihps EM430ST TEM used is housed in a triple- 
bi/story building designed for low vibration, and provides 
continuous contrast transfer to ^ 1/(0. 19nm) at Scherzer 
defocus. It is equipped with a ±15° side-entry goniome- 
ter specimen stage. A Gatan double tilt holder enables 
±10° tilt around an orthogonal tilt axis. The largest 
orientation difference which can be achieved using this 
double tilt holder in the microscope is therefore 35.6° 




B. Determining the angle of effective tilt projected 
onto an image 

In order to establish the spatial relationship between 
reciprocal lattice vectors inferred from images taken at 
different specimen tilts, the direction of the tilt axis rel- 
ative to those images must be known. The tilt axis di- 
rection is defined via the right hand rule, as orthogonal 



to the relative motion of parts of the specimen as the 
goniometer reading is increased. In a single tilt, the axis 
is perpendicular to the electron beam and parallel to the 
micrographs. This is true also of the effective tilt axis 
in a double tilt, provided the two specimen orientations 
are symmetric about the zero tilt position. We limit our 
discussion of double tilts to this case. 

We determined the projection of both tilt axes of a 
Gatan double tilt holder onto TOOK HREM images by 
examining Kikuchi line shifts during tilt in the 1200mm 
diffraction pattern of single crystal silicon, and then cor- 
recting for the rotation between that diffraction pattern 
and the image To be specific, with a micrograph 
placed in front of the operator with emulsion side up as 
in the microscope, with zero azimuth defined as a vector 
from left to right, and with counterclockwise defined as 
the direction of increasing azimuth, the projection of the 
Ti axis on 1200mm camera- length diffraction patterns in 
our microscope is along 114.0°. The rotation angle be- 
tween electron diffraction patterns at the camera length 
of 1200mm and TOOK HREM images is 42.9°. Therefore 
the direction of the projection of Ti on TOOK HREM im- 
ages is along —156.9°. The direction of the projection of 
the second tilt axis, T2, on TOOK HREM images is along 
113.1°. 



C. A reference coordinate system 

We then consider a coordinate system fixed to the mi- 
croscope, for measuring reciprocal lattice vectors from 
the power spectra of TOOK HREM images. The y and z 
directions are defined to be along — Ti and the electron 
beam direction, respectively, as shown in Figure ^. The 
projection of these tilts on the power spectrum of a TOOK 
HREM image is shown in the 2nd inset of Figure ||. Az- 
imuths in the remainder of this paper are all measured in 
the xy plane of this coordinate system, with the x or T2 
direction set to zero. Because the Ti direction is defined 
in our coordinate system as the negative y-direction, az- 
imuthal angles are measured on micrographs from a di- 
rection 90° clockwise from the Ti direction, when the 
emulsion side is up. 



D. Double tilting 

The specimen was first tilted about T2 to z92 — 9.T° 
while 1^1 remained at 0°, made eucentric, then tilted 
about Ti to di = 15.0°. The first HREM image was taken 
at this specimen orientation of (i9i = 15°,i?2 = 9.T°). A 
similar sequence was applied to take a second HREM im- 
age at the second specimen orientation of (i?i — —15°, 
'!?2 = — 9.T°). The process can be modeled with a simple 
matrix calculation . 

Because of the importance of repeatable quantitative 
tilts, effects of "mechanical hysterisis" were minimized 
by inferring all relative changes in tilt from goniometer 
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Projection of Tilt Axes 
On Electron Micrograph 



FIG. 2: A schematic illustration of the coordinate system 
set-up for measuring reciprocal lattice vectors. The coordi- 
nate system is fixed to the microscope. The y and z axes 
are defined to be along -Ti and the electron beam direction, 
respectively. The projections of Ti and T2 on 700kx magni- 
fication micrographs, as well as of Ti on 1200mm diffraction 
patterns, are shown as well. 



TABLE I: The g-spacings [nm^^] and azimuths [°] of three 
spots measured from the power spectra of images at two tilts, 
as well as calculated Cartesian coordinates [nm~^] of those 
reciprocal lattice spots in a common coordinate system. This 
constitutes a minimal data set for analyzing the lattice in 
three- dimensions. 

Spot n Qm "dl "02 Qx Qy flz 

1 4.73 79.2 -15.0 -9.7 0.861 4.53 -1.01 

2 4.77 -11.6 -15.0 -9.7 4.52 -1.15 -1.03 

3 4.14 32.6 4-15.0 -f9.7 3.37 2.04 1.27 



of 330°C ||T5[. The specimen was disk-cut, abraded from 
the glass substrate side and dimpled by a Gatan Model 
601 Disk Cutter, a South Bay Technology Model 900 
Grinder and a Gatan Model 656 Precision Dimpler, re- 
spectively. The specimen was then argon ion-milled by a 
Gatan DuoMill for about 5 hours to perforation prior to 
the TEM study, at an incidence angle of 3°. 



IV. EXPERIMENTAL RESULTS 

A. Diffraction from a known to check column 
geometry 

Calibration of these algorithms with geometry in our 
microscope was first done with diffraction data from a Si 
crystal. Diffraction patterns of (100) Si along the [116] 
and [116] zone axes were obtained by tilting about Ti and 
T2. The lattice parameters determined are (a = 0.383nm, 
b = 0.387nm, c = 0.386nm, a = 60.0°, /3 = 119.6°, 
7 = 119.1°). This set of chosen basis defines the rhombo- 
hedral primitive cell of the Si f.c.c lattice. Compared with 
the literature values of Si (a = 0.384nm, b = 0.384nm, 
c = 0.384nm, a = 60°, /3 = 120°, 7 = 120°), the angular 
disagreements are less than 1° and spatial disagreements 
are less than 1%. These uncertainties are comparable to 
those obtained by other techniques of submicron crystal 
analysis P, |l|, jll, IH, [l7 . 



readings taken with a common direction of goniometer 
rotation. The rotation sequences were "initialized" by 
first tilting past the starting line, and then returning to it 
in the direction of subsequent motion. Nonetheless both 
the precision of angle measurement, and our inability to 
observe lattice fringes during rotation, were shortcom- 
ings that microscopes designed to apply these strategies 
routinely must address. 



E. Specimen preparation 

The tungsten carbide thin film was deposited by 
PECVD on glass substrates by introducing a gaseous 
mixture of tungsten hexacarbonyl and hydrogen into a 
RF-induced plasma reactor at a substrate temperature 



B. Nanocrystal images to infer lattice parameters 
of an unknown 

The micrographs in Figure ^ show a nanocrystal in 
a film rich in tungsten carbide, at the orientations of 
= 15°, ^2 = 9.7°) and {di = -15°, i?2 = ~9.7°) 
respectively. The coordinates of periodicities in micro- 
graph power spectra, as well as in the common reference 
coordinate system, are listed in Table |, in much the same 
format as is diffraction data for stereo analysis [^ . 

1. Matching the lattice: 

The lattice spacings and inter-spot angles of periodic- 
ities in image power spectra were used to look for con- 
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FIG. 3: The tilt protocol for inferring the lattice structure of f.c.c. WCi-x by viewing a WCi-x crystal from its [001] zone 
and [112] zone, along which the two most easily resolved lattice plane sets, the {200} and {111} lattice planes, will show lattice 
fringes in HRTEM images. The two zones are 35.3° apart, within the tilting limit of the microscope with a Gatan double tilt 

holder. The crystal is to be viewed along its [001] zone in the first specimen orientation, then tilted along the (220) lattice 
planes to the [112] zone. The (220) lattice planes must therefore be perpendicular to the effective tilt axis. The projections 
along the [001] and [112] zones together with the effective tilt axis have been drawn so that their azimuths are consistent with 
those in the 700K HRTEM images. The experimental HRTEM images of a WCi-a; nano-crystal acquired using this protocol 
and their power spectra are shown in the bottom. The actual tilt sequence is to tilt along -T2 by 19.5° followed by tilting along 
-Ti by 30.0°, where Ti and T2 denote the side-entry goniometer tilt axis and the second tilt axis of a Gatan double tilt holder, 
respectively. 
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TABLE II: The results of a three-dimensional match of mea- 
sured spacings [nm] and interspot angles [°], with those pre- 
dicted from the literature (denoted with a caret) for the face- 
centered cubic crystal WCi-^,. This phase is the only one 
from a list of 36 tungsten carbides and tungsten oxides whose 
predicted spacings and interspot angles agreed the measure- 
ments within a tolerance of 1.5% and 1.5°, respectively. 

Spot n d„ ipij^n (hkl) dhki ^[%] 'fiiji^n Sip 



0.212 54.2 (200) 0.212 0.5 54.7 0.5 
0.209 56.2 (020) 0.212 1.4 54.7 1.5 
0.242 90.8 (llT) 0.245 1.2 90.0 0.8 



sistent indexing alternatives from a set of 36 tungsten 
carbide and oxide candidate lattices including WCi-x- 
When an angular tolerance of 2° and a spatial tolerance 
of 2% are imposed, only WCi^x provides a consistent 
indexing alternative. As summarized in Table ||, the 
Miller indices of the three observed spots then become 
(200), (020) and (llT). The images of Figure | thus rep- 
resent WCi_2, [001] and [112] zones, respectively. The 
azimuth of the reciprocal lattice vector (220), 



'^(220) 



("^(200) + [180° + </'(020)]) 



123.. 



(7) 



deviates from the projection of the effective tilt axis by 
only 0.3°. Therefore the (220) lattice planes are perpen- 
dicular to the effective tilt axis as per Figure ||, and the 
data acquired are consistent with the expectation for fee 
crystals outlined in Figure The actual tilting path in 
the Kikuchi map of crystal A is shown there as well. 

From the indexing suggested by this match, the a;, y, 
and z coordinates of the reciprocal lattice basis vectors 
a*, b\ and c* may be inferred. This is shown in Ta- 
ble m, along with the resulting lattice parameters and 
a comparison with literature values. The resulting er- 
rors in a and b are less than 1.3%, while the error in c 
(which is orthogonal to the plane of the first image) is 
larger (around 2.3%). Both because of tilt uncertainties 
and reciprocal lattice broadening in the beam direction, 
uncertainties orthogonal to the plane of the specimen are 
expected to be larger than in-plane errors S. 



2. Building a triplet from scratch: 

By generating linear integral combinations of the mea- 
sured periodicities in reciprocal space (i.e. vector triplets 
of the form nigi + 712.92 + ^^333 where the rii are integers) 
until a minimal volume unit cell is obtained (there will 
be more than one way to achieve the minimum) , a primi- 
tive triplet for the measured lattice can be inferred quite 
independent of any knowledge of candidate lattices. The 
primitive cell parameters determined are also listed in 
Table [II. With respect to literature values, these show 



spatial disagreements of less than 1.6%, and angular dis- 
agreements of less than 1.5°. Although inference of the 
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FIG. 4: Illustration of all ways (from Figure 0) to verify the 
three dimensional lattice parameters of fee and bcc crystals 
from a pair of lattice images, given an ability to image only 
lattice spacings down to half the unit cell side, and a tilt range 
of less than 60 degrees. The image pairs are rotated so as to 
illustrate the direction of tilt between the two zones (crystal 
orientations) involved. The protocols shown with double ar- 
rows provide access to cross-fringes (and hence information 
on which direction to tilt) at both ends of the tilt sequence. 
Also note that the first and last entries in each column sim- 
ply provide low and high resolution views, respectively, of the 
same experiment. 



TABLE III: Primitive (top) and face-centered (bottom) recip- 
rocal lattice basis triplets (spacings in nm, angles in degrees) 
inferred from our two HREM images of crystal A, along with 
a comparison of the lattice parameters for each cell which 
follow therefrom. The indexing of the observed reciprocal 
lattice vectors in these two cases was inferred directly for the 
primitive cell, by minimizing cell volume, and by matching 
to "textbook" parameters for WCi_a: in the case of the face- 
centered cell. Except for the choice of basis triplet, the two 
measured cells refer to exactly the same inferred lattice. 

measured 0.298 0.299 0.296 120.0 58.7 119.8 
"book" 0.300 0.300 0.300 120 60 120 



params 
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"book" 0.425 0.425 0.425 
measured 0.424 0.419 0.415 



90 90 90 
88.5 90.8 89.3 



"conventional cell" from the primitive cell alone is possi- 
ble, the process has not been attempted here because of 
complications attendant to measurement error. 
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3. Phase Identification 

Determining a reciprocal lattice triplet, and inferring 
lattice parameters therefrom, are of course not equivalent 
to confirming the existence of a particular phase. In or- 
der to draw a more robust conclusion about the makeup 
of crystal A, we extended our analysis of the structure 
to other lattices capable of indexing the observed spots, 
albeit with larger errors in spacing and interspot angle. 
When the spatial and angular tolerances of our candi- 
date match analysis are increased to 3° and 3%, there 
are many tungsten oxide and carbide candidates in ad- 
dition to WCi-x which show consistent lattice spacings 
and inter-planar angles [ p^ . 

In order to eliminate these candidates, it was necessary 
to confirm, using power spectra of amorphous regions in 
each image, that the spatial frequencies in Figure ^mages 
were continuously transferred within the first passband 
lH, By then assuming that projected reciprocal lat- 
tice frequencies make their way into the exit surface wave- 
field (at least at the thin edges of the particle), all the 
candidates except WCi^x are eliminated. Specifically, it 
was found that for each of the candidates except WCi_a;, 
along one or more of the suggested zone axes at least one 
reciprocal lattice vector shorter than the experimental 
one(s) is missing in a power spectrum p^ . An example 
of this is the match with hexagonal WC^, (a = 1.058nm, 
c = 1.335nm). In this case the Miller indices suggested 
for spot 3 (422) were inconsistent with the fact that the 
(211) is absent from the power spectrum of the image on 
the right side of Figure ^. 

Our conclusion that "this crystal is WCi-a;" (and as 
we see later most of the other crystals in the specimen) 
is consistent with knowledge of the formation conditions, 
as well as with X-ray powder and EDS analysis of the 
same film. 



C. The effective tilt direction, and recurring fringes 

In addition to serving as a guide for correctly choosing 
the azimuth of the crystal before tilting between desired 
zones, knowledge of the tilt axis direction plays another 
role: that of highlighting lattice fringes present in both 
specimen orientations, but caused by one and the same 
set of lattice planes. 

In single tilt experiments, the tilt axis is simply Ti. 
This is always perpendicular to the electron beam and 
hence parallel to the micrographs. Any reciprocal lattice 
vector parallel or antiparallel to Ti remains in Bragg con- 
dition throughout the whole tilting process, regardless of 
the amount of tilt t9i. If the spacing is large enough to 
be recorded in the images, the same lattice fringes are 
seen perpendicular to the projection of Ti in any HREM 
image as well. 

For double tilt experiments, it is convenient to intro- 
duce the concept of an effective tilt axis. The effective 
tilt axis is analogous to the tilt axis in a single tilt ex- 



periment, in that the double tilt can be characterized by 
a single tilt around the effective tilt axis of angular size 
equal to that in the double tilt. This effective axis is 
perpendicular to the electron beam and hence parallel to 
the micrographs only if the two specimen orientations are 
symmetric about the untilted position. 

Considering only double tilts falling into this category, 
let (i?i, 1^2) and (— —'02) denote the specimen orienta- 
tions before and after. The effective tilt axis direction has 
an azimuth (with respect to our reference x-direction) of 



(fieff = tan" 



sin(i9i) 
tan(i92) 



(8) 



A proof of equation ^ is given in Appendix A. There 
exists a 180° ambiguity in the direction of the effective 
tilt axis using equation]^. This ambiguity can be resolved 
through the knowledge of the actual tilting sequence. In 
our experiment -di = 15°, i?2 = 9.7°, (peff = 123.5°. This 
is the effective tilt axis direction mentioned in previous 
sections. 

Lattice planes perpendicular to the effective tilt axis, 
in the double tilt case, diffract and are visible at initial 
and final, but not intermediate, specimen orientations. 
This result inspired further experimental work on, and 
modeling of, fringe visibility loss during tilt [Q. One re- 
sult of this exercise was a prediction that 0.213nm fringes 
deviating by as much as 4° from the effective tilt direc- 
tion in a lOnm WCi-x specimen will remain visible after 
a 35.6° tilt. This was confirmed by experiment on these 
specimens 1). The result in turn serves to constrain the 
probability and error analyses below. 



D. Tilt limitations and chances for success 

This section addresses the chances for successful 3D 
cell determination from images, depending on properties 
of both microscope and specimen. Such matters are con- 
sidered in more detail elsewhere 

In a microscope with a single-axis tilt of at least ±35.3° 
and a stage capable also of 180° rotation, any cubic crys- 
tal with an [001] zone in the beam direction at zero tilt 
can be re-aligned by azimuthal rotation until its (220) re- 
ciprocal lattice vector is parallel to Ti. With an untilted 
tilt-rotate stage, this would allow a crystal's [001] zone to 
remain aligned with the beam throughout the rotation. 
Subsequent tilting by 35.3° will lead to the [112] zone, 
and the lattice structure in three dimensions confirmed 
ala Figure |[ 

Under these conditions, any cubic crystal showing [001] 
zone cross fringes can be tilted so as to reveal a third 
spacing. Hence the probability of success with any given 
crystal is that of finding a randomly-oriented crystal ori- 
ented with [001]-zone fringes visible. Fortunately for this 
method, the spreading of reciprocal lattice points due to 
finite crystal thickness t allows one to visualize fringes 
within a half angle 0( of order 1/t surrounding the exact 
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Bragg condition. Otherwise, cross-fringes would be rare 
indeed! 

The sohd angle subtended by this visibility range for 
lattice planes intersecting along the [001]-zone allows us 
to calculate the probability that a randomly-oriented 
crystal will show cross-fringes of specified type. For ex- 
ample if we approximate the cross-fringe region with a 
conical bundle of directions about each zone, then for the 
special case of spherical particles the fraction of crystals 
showing the fringes of zone x is: 



arccos 



9d- gf -'2gx9i, 



Px — n{l — cos[Qt]), where Ot — arcsin[ 



9t -9l + '^9x9t ^ 
V^9x9\ 



n is multiplicity of zone x (e.g. n = 'i for cubic x =[001]), 
gx = l/d, gx = 1/A, and gt = f /t, where t is the thick- 
ness of the crystal in the direction of the beam and / is 
a parameter of order one that empirically accounts for 
signal-to-noise in the method used to "visualize" fringes. 
For example, we expect / to decrease if an amorphous 
film is superimposed on the crystals being imaged. The 
half-angle Qt is a pivotal quantity in both the probability 
and accuracy of fringe measurement. 

A plot of the probability for seeing [001] cross- fringes of 
spacing d = 0.202nm, as a function of specimen thickness 
t for both spherical and laterally-infinite (rel-rod) parti- 
cles, is shown in Fig. ^. Here we've used fit parameter 
/ = 0.79 based on data points (also plotted) that were 
obtained experimentally for particles of varying size from 
HREM images of Au/Pd evaporated onto a carbon film. 
As you can see, the probability of encountering cross- 
fringes improves greatly as crystallite size decreases to- 
ward a nanometer. Of course, as discussed in the next 
section, this "reciprocal lattice broadening" is accompa- 
nied by a decrease in the precision of measurements for 
individual lattices. 

Due to the tilt limits of the specimen holder in our mi- 
croscope, the first HREM image along the [001] zone of 
a WCi-x nano-crystal had to be taken at a nonzero i^i 
orientation. Azimuthal symmetry is thus broken. Our 
solution was to find a [001] nano-crystal whose (220) re- 
ciprocal lattice vector was by chance parallel to the ef- 
fective tilt axis, then tilting to the 2nd orientation. Thus 
nano-crystal A was identified (by coincidence) to have an 
appropriate azimuth during real time study of the (200) 
and (020) lattice fringes. Tilting by 35.3° was done there- 
after. 

For the probability of success in our case, we must mul- 
tiply Px by the probability of viewing (111) fringes after 
tilting a [001] crystal with random aziumth by 35.26°. 
This probability of finding a 3rd spacing takes the form 
P3 — m(5/7r, where m is the multiplicity of target spac- 
ings (e.g. m = 4 for a four-fold symmetric [001] starting 
zone), and the "azimuthal tolerance half-angle" 5 (again 
in the spherical particle case) obeys the implicit relation: 

^ .tan^o cos 7 , 

Oo = arctan + arctan^-g — — — , 

cos sm 7 - (cos Oo sm Sy 

where 9o is the required tilt (in our case 35.26°) and 



The probability is also plotted as a function of speci- 
men thickness in Fig. |^, along with the product PxPs- 

These models predict a probability of success with the 
strategy adopted in our experiment, for the "large" lOnm 
WCi-a; crystals in our specimen, of pxPs = 7.3 x 10^* x 
0.371 = 2.4 X 10"^. Hence only one in every 1500 crystals 
will show [001] cross fringes, and one in every 4000 will 
be suitably oriented for 3D lattice parameter determina- 
tion. This is consistent with our experience: The image 
of crystal A was recorded in one negative out of 22, each 
of which provided an unobstructed view of approximately 
100 crystals. 

As mentioned above, using a microscope capable of 
side-entry goniometer tilting by ±35.3° with a tilt-rotate 
stage, the 3D parameters of all cubic crystals, when un- 
tilted showing [001] zone cross-fringes, could have been 
determined. According to Fig. 14a, the fraction of parti- 
cles 2 nm in thickness that are oriented suitably for such 
analysis approaches 1 in 100. Moreover, with a goniome- 
ter capable of tilting by ±45° plus computer support for 
automated tilt/rotation from any starting point, each 
unobstructed nano-crystal in the specimen could have 
been subjected to this same analysis after a trial-and- 
error search for accessible [001] zones. Thus a significant 
fraction of crystals in a specimen become accessible to 
these techniques, with either a sufficient range of precise 
computer-supported tilts, or if the crystals are sufficiently 
thin. 

Subsequent work [ |20| has shown that information on 
tilt protocols and fringe visibility for crystals of a given 
thickness can be elegantly sumarized with spherical maps 
like those in Figure |. These are a direct-space analog to 
reciprocal-space Kikuchi maps, in which band thickness 
is proportional to d (rather than 1/d). If the sphere used 
has a diameter equal to specimen thickness, then the first 
order effect of changing thickness simply increases the 
separation between zones while holding the width of the 
bands fixed. 



V. PITFALLS AND UNCERTAINTIES 



A. Cautions involving specimens and contrast 
transfer 



In this section, we discuss effects warranting caution. 
In the next section, models of lattice parameter uncer- 
tainty are discussed. 

High electron beam intensities can cause lattice rear- 
rangement in sufficiently small nanocrystals, as well as 
changes in the orientation of a thin film (e.g. due to 
differential expansion). Sequential images of the same 
region at fixed tilt might allow one to check for such 
specimen alterations. 
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FIG. 5: A plot of the fraction of randomly oriented grains of WCi^x showing [001]-zone cross fringes (px), and the fraction of 
such cross-fringe grains oriented so that a random tilt of 35.3 ° will allow imaging of a 111 periodicity (ps), as a function of 
specimen thickness in the direction of the electron beam. IiqI 



Loss of periodicities in the recorded image, due to 
damping and spherical aberration zeros, were discussed 
in the section above on experimental design. Nonethe- 
less, careful observations of more than one crystal, and 
image simulation as well, may be useful adjuncts when- 
ever this technique is applied. We illustrate this below, 
with a "two-dimensional" experiment done to assess the 
size of errors due to finite crystal size and random orien- 
tation. The result is of help in the section on modeling 
uncertainties that follows. 

A recent paper on HREM image simulations |l^ indi- 
cated that deviations in orientation of a 2.8nm spherical 
palladium nano-crystal from the zone axes may result in 
fringes unrelated to the structure of the particle. Vari- 
ability in measured lattice spacings was also reported to 
be as high as several percent, with the highest reach- 
ing 10%. To compare such results with our experimental 
data, 23 single crystals free of overlap with other crystals. 



and each showing cross-fringes, were examined. The pro- 
jected sizes of these crystals range from 3.7nmx3.8nm to 
lO.Snmx 7.8nm. The spacings and angles between fringes 
are plotted in Fig. |^. 

Observed cross-fringes in the HREM images fall into 
two categories, according to their spacings and angles. 
The first category is characterized by a 90° inter-planar 
angle between two 2.12A lattice spacings. The second one 
by two inter-planar angles of 55°, 70° and three lattice 
spacings of 2.12A, 2.121, 2.441. Only the spacings of 
2.441 and 2.121 and the corresponding angle of 55° have 
been shown in Figure ^. Two conclusions can be drawn. 

First, since the two categories of cross- fringes match 
those along the [001] and [Oil] zone axes of WCi_a;, the 
only two zones which show cross-lattice fringes in our 
HREM images, the thin film consists mainly of WCi^j,. 
X-ray powder diffraction work on the film supports this 
conclusion p^ . 



11 




FIG. 6: Visibility maps for the two largest fringes visible from body-centered, face-centered, and diamond face-centered lattices. 
Sphere diameter (relative to band width and fringe spacing) is proportional to specimen thickness, here chosen to be about 
5 times the cubic cell side a. Note the dominance of crossed (110) fringes at the three-fold (111) zone in the body-centered 
case, the dominant crossed (111) fringes at the two-fold (110) zone in the face-centered cases, and the wider disparity between 
largest and next-to-largest spacings when the diamond glide is added to the lattice. 
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FIG. 7: The spacings and interplanar angles measured from the cross lattice fringes of 23 nano-crystals free of overlap with 
other crystals in the HRTEM images. Insets (b) and (c) are magnified plots of the two regions in (a) where all the data are 
concentrated. The specific combinations of lattice spacings and interplanar angles corresponds to the [001] and [Oil] zone 
images of WCi-a; and hence indicate WCi-a, being the only present phcise in the thin film. The meEisured lattice spacings have 
a standard deviation from their mean of less than 1.5%, and interplanar angles have a standard deviation of 1.2°. 
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Secondly, for nano-crystals free of overlap with other 
crystals, the observed lattice spacings and interplanar an- 
gles in HREM images have a standard deviation from 
the mean of less than 1.5%, and a standard deviation of 
less than 1.3°, respectively. We have not observed any 
seriously shortened or bent fringes. Nonetheless we rec- 
ommend that such fringe abundance analyses go hand 
in hand with stereo lattice studies of nano-crystal speci- 
mens, and that comparative image simulation studies be 
done where possible as well. 



B. Uncertainty forecasts 

Earlier estimates as well as the typical size of 
diffraction broadening in the TEM, suggest that lattice 
parameter spacing errors may in favorable circumstances 
be on the order 1%, and angle errors on the order of 
1°. Experiment, and a more detailed look at the theory 
, now support this impression. We will focus the dis- 
cussion here on equant or spherical nanocrystals. The 
results should be correct within 10% for other (e.g. thin- 
foil) geometries of the same thickness. 

The three sources contributing to the lattice spacing 
measurement uncertainty in images are: expansion of the 
reciprocal lattice spot in the image plane, uncertainty in 
the camera constant, and expansion of the reciprocal lat- 
tice spot along the electron beam direction, in order of 
decreasing relative effect ||^ . The uncertainty in our cam- 
era constant is measured to be about 0.5%. Uncertainties 
from the first and third sources above are on the order of 
1% and 0.01%, respectively, for a typical lattice spacing 
of 0.2 nm. The in-plane/out-of-plane error ratio is on the 
order of 10. 

Sources contributing to the measurement uncertainty 
of lattice parameters along the electron beam direction, 
when the specimen is un-tilted, include uncertainty in go- 
niometer tilt as well as sources analogous to those above. 
Observation of reciprocal lattice vectors further out of the 
specimen plane (i.e. of fringes at high tilt) reduces the 
measurement uncertainty of "out-of-plane" parameters. 
The measurement uncertainty of inter-planar angles in 
images is due to lateral uncertainty in their associated 
reciprocal lattice spots in the image plane. 

Using a mathematical model of these errors |^ , we pre- 
dict spacing uncertainties in a 10 nm nanocrystal, tilted 
by ±18°, of 2.1% for an imaged spacing and of 8.6% for 
a lattice parameter perpendicular to the plane of the un- 
tilted specimen. This large "out-of-plane" uncertainty is 
a result of the small tilt range available with our high 
resolution pole piece. The estimated interplanar angle 
uncertainty is about 2.3°. These predicted uncertainties 
are between 2 and 3 times the errors observed here, 
and hence of the right order of magnitude. 

The model suggests that lattice parameter uncertain- 
ties will decrease as camera constant and tilt uncertain- 
ties decrease, and will also decrease as the tilt range used 
for the measurement increases. It suggests that the lat- 



tice parameter errors will increase as crystal thickness 
goes down. The ease of locating spacings, however, goes 
up as crystal thickness decreases. Hence the best candi- 
dates for application of the protocols here may be crystals 
in the 1 to 20 nm range. Improved tilt accuracy (hope- 
fully with computer guidance), and low- vibration tilting 
so that fringes may be detected as orientation changes, 
would make these strategies more accurate and widely 
applicable as well. 



VI. CONCLUSIONS 

When considered from the perspective of direct space 
imaging, crystals offer a discrete set of opportunities for 
measuring their lattice parameters in three dimensions. 
Ennumerating those opportunities for candidate lattices, 
or lattice classes, opens doors to the direct experimental 
determination of nano-crystal lattice parameters in 3D. 
A method for doing this, and lists of those opportunities 
for the special case of cubic crystals, are presented here. 

We apply this insight to inferring the 3D lattice of a 
single crystal from electron phase or Z-contrast images 
taken at two different orientations. For nano-crystals in 
particular, a double-axis tilt range of ±18° degrees al- 
lows one to get such data from all correctly-oriented cubic 
crystals with appropriate spacings resolvable in a pair of 
images taken from directions separated by 35.3°. In the 
experimental example presented, we find less than 1.5% 
spatial and 1.6° angular disagreements between the in- 
ferred primitive cell lattice parameters of a 10 nm WCi-x 
nano-crystal, and literature values. 

We further present data on the variability of lattice 
fringe spacings measured from images of such randomly- 
oriented lOnm WCi_2; crystals in electron phase contrast 
images. The results suggest that measurement accura- 
cies of 2% in spacing and 2° in angle may be attainable 
routinely from particles in this size range. Smaller size 
crystals may be easier to obtain data from, but show 
larger uncertainties, while larger or non-randomly ori- 
ented crystals (especially if guesses as to their structure 
are unavailable) may be more challenging to characterize 
in three dimensions. 

Precise knowledge of the tilt axes, as projected on the 
plane of a micrograph, is crucial to implementation. This 
information, if coupled with on-line guidance on how to 
tilt from an arbitrary two-axis goniometer orientation in 
any desired direction with respect to the plane of an im- 
age or diffraction pattern, could make this strategy and 
related diffraction strategies for lattice parameter 
measurement more routine. Future microscopists might 
then be able to interface to individual nanocrystals much 
as the nano-geologist in the first paragraph of this paper 
examined her "hand specimen" . 

Lastly, diffraction can also be used in this "stereo anal- 
ysis" mode (with crystals large enough to provide diffrac- 
tion patterns), although the easier accessibility of high 
spatial frequencies via diffraction often makes the large 
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tilts used here unnecessary. They can also be put to 
use in darkfield imaging applications, by forming images 
of the specimen using "beams" diffracted by the period- 
icities which serve as diagnostic of a given lattice (for 
example those associated with each of the protocols of 
Figure [|) . 

Images so taken of nano-crystalline specimens, for ex- 
ample at two tilts with three different darkfield condi- 
tions, would be expected to show correlations among 
that subset of the crystals correctly oriented for diffrac- 
tion with all three reflections. Although this strategy 
may never allow precise lattice parameter determinations 
given limits on objective aperature angular size, it may 
be a very efficient way to search for crystals correctly- 
oriented and of correct type for one of the imaging pro- 
tocols described here. Moreover, because such lattice- 
correlations in three dimensions contain information be- 
yond the pair-correlation function, they may be able 
to support the new technique of fluctuation microscopy 
pl| , p2 , [23| in the study of paracrystalline specimens like 
evaporated silicon and germanium [ p4| , whose order- 
range is too small for detection by other techniques. 



APPENDIX A: EFFECTIVE TILT AXIS 
AZIMUTH 

Let {-di, 1^2) and {"S'l, '02) denote orthogonal tilt val- 
ues for two specimen orientations, and 4'eff the azimuth 
of the effective tilt axis between these orientations. Any 
reciprocal lattice vector with untilted Cartesian coordi- 
nates I (7), and with identical micrograph coordinates \gm) 
at the two tilted orientations, will (following equation H) 
obey: 



where 



/ g cos (ipeff) 
3 sin {(fieff) 



(A2) 



V 







Expanding, this gives three equations which, combined 

be solved for 



with equation A2 



the three unknowns i?']^, §'2, ^effi to get: 



^'1 = -^1, 



-'dn, and 



LPeff = tan" 



sin(i9i) 
tan(i92) 



(A3) 



(A4) 



^(^i,i?2)|g™> = 1.9) -^(^'i,^9'2)|.9. 



(Al) 



This provides an equation for the azimuth of the effec- 
tive tilt, and confirms that symmetry about the zero tilt 
position is a necessary and sufficient condition for the 
reciprocal lattice vector jg), and it's associated lattice 
fringe, to show a common direction in micrographs at 
both tilts. 
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